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ABSTRACT 

We present detailed analysis of three globular cluster X-ray sources in the XMM-Newton extended 
survey of M31. The X-ray counterpart to the M31 globular cluster Bo 45 (XBo45) was observed 
with XMM-Newton on 2006 December 26. Its combined pn+MOS 0.3-10 keV lightcurve was seen to 
vary by ~10%, and its 0.3-7.0 keV emission spectrum was well described by an absorbed power law 
with photon index 1.44±0.12. Its variability and emission is characteristic of low mass X-ray binaries 
(LMXBs) in the low-hard state, whether the accretor is a neutron star or black hole. Such behaviour 
is typically observed at luminosities <10% Eddington. However, XBo45 exhibited this behaviour at 
an unabsorbed, 0.3-10 keV luminosity of 2.5±0.2xl0 38 erg s" 1 , or ~140% Eddington for a 1.4 M Q 
neutron star accreting hydrogen. Hence, we identify XBo45 as a new candidate black hole LMXB. 
XBo45 appears to have been consistently bright for ^30 years, consistent with theoretical prediction 
for a globular cluster black hole binary formed via tidal capture. Bo 375 was observed in the 2007, 
January 2 XMM-Newton observation, and has a two-component spectrum that is typical for a bright 
neutron star LMXB. Bo 135 was observed in the same field as Bo 45, and could contain either a black 
hole or neutron star. 
Subject headings: X-rays: binaries — galaxies: individual: M31 — methods: data analysis 



1. INTRODUCTION 

Of the 14 bright low mass X-ray binaries (LMXBs) 
associated with Galactic globular clusters (GCs), 13 
are known to contain neutron stars, and the other one 
probably does as we ll (see e.g. IWhite fc Angelim] 120011 : 
lin't Zand et alJl2004L and references within). However, 
black hole binaries have been seen in extragalactic GCs, 
provin g that they form i n suc h environments. For ex- 
ample, [MaccaroneleF^i] (|2007h discovered the first GC 
black hole binary in the giant elliptical galaxy NGC 4472. 

The lack of known GC black hole binaries has been 
a subject of g r eat in terest for theoretical modellers. 
iKalogera et al.l (|2004f ) show that black hole binaries 
formed through exchange interactions should have duty 
cycles of ^0.001, consistent with the absence of GC black 
hole binaries at the time. They also predict that black 
hole binaries formed through tidal capture of a main se- 
quence star would be bright, persistent X-ray sources, 
and infer from their absence that tidal capture proba- 
bly disrupts the main sequence star. However, neither 
of these results preclude the discovery of new GC black 
hole binaries. 

The bright X-ray sources of M31 have been studied 
for o ver 25 yea rs, with Einstein 
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Shaw Greening et all 120081) . 



iDi Stefano et al.l (|2002D conducted a Chandra survey 
of selected regions of M31, and found that most of their 
bright X-ray sources were associated with GCs, with 

10% of GC sources exhibiting 0.5-7 keV luminosities 
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They also showed the M31 GC popula- 



e.g. IDi Stefano et alj|2002t lKaaretJl2002t IWilliams et al 
I2004T ) and XMM-Newton (see e.g. IShirev et al" 
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tion to be significantly different from that of the Milky 
Way as ~30% of the M31 GC X-ray sources exhib- 
ited luminosities >10 37 erg s _1 , whereas only one out 
of 11 Galactic GC X -ray sources exceeded 10 37 erg s _1 
(IVerbunt et alj|1995l) . These results contradict those of 
ISupper et al.l (|200lf) . who concluded that the GC X-ray 
p opulations of the Milky W ay and M31 were similar. 

Di Stefa no et al.l (|2002f ) also compared the optical 
properties of M31 GCs with and without X-ray sources, 
as well as M31 GCs and Milky Way GCs with X-ray 
sources. They compared optical colours and apparent 
magnitudes, radial velocities, metallicities and colour ex- 
cesses, as well as core radii. They found the only possi- 
bly significant difference to be in their luminosities: the 
median luminosity of M31 GCs with X-ray sources was 
^0.55 magnitudes higher than that of M31 GCs with- 
out. Similarly, M31 X-ray GCs were slightly brighter in 
V than Milky Way X-ray GCs. IDi Stefano eTaLl (pfM ) 
speculated that this might associate X-ray sources with 
higher GC masses. 

We report on three bright X-ray sources in M31, all 
associated with populous, old globular clusters: Bo 45, 
Bo 135 and Bo 375. We will argue that Bo 45 contains 
a good black hole (BH) candidate, while Bo 375 is likely 
to contain a neutron star (NS). Bo 135 is more likely to 
contain a NS, but we cannot rule out a BH accretor. 

We first review the emission and variability of LMXBs 
in their various states in Section O and then discuss 
known properties of our targets and their host clusters 
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in Section [3j We describe the observation and analysis 
in Section 01 and our results in Section [5j We finally 
present our discussion and conclusions in Section [5] 

2. VARIABILITY AND EMISSION SPECTRA OF LMXBS IN 
DIFFERENT STATES 



It was shown bv lvan der Klisl (J1994D that NS and BH 
LMXBs accreting at low rates exhibit remarkably sim- 
ilar emission spectra and variability. In this low-hard 
state, the power density spectrum (PDS) may be char- 
acterised by a broken power law with spectral index 7 
changing from ~0 to ~1 at some break fre quency; also, 
the r . m.s. variability is high (~10-50%, Ivan der Klisl 
fl99J [19951 iMcClintock fc Remillardl l200l . However, 
at higher accretion rates, the PDS may be described 
by a simple power law with 7 ~1 and the r.m.s . vari- 
ability is only a few percent (|van der Klisl I1994L 119951 : 
IMcClintock fc Remillardl l2006h . 

This change in variability from low to high accretion 
rate is accompanied by a change in the emission spec- 
tra of LMXBs: low-state LMXBs have hard spectra 
that are characterised by power laws with photon index 
—1.4-1.7, regardless of the pri mary (jvan der Klisl 119941 
IMcClintock fc Remillardl I2006T ) . Higher accretion rates 
lead to distinctive emission spectra, depending on the 
accretor. BH LMXBs exhibit a thermally dominated 
state (also known as the high/soft state) where a 0.7- 
2 keV dis c blackbody contributes >90% of the 0.3-10 
keV flux (|McClintock fc Remillardl 12006): they also ex- 
hibit a "steep power law" state, where an additional 
power law component is observed, with a photon in- 
dex > 2.4 (|McClintock fc Remillardl 120061) . Meanwhile, 
the emission spectra of high accretion rate NS LMXBs 
remain predominantly non-thermal, although a second, 
blackbody component becomes increasingly important at 
higher l uminosities, contributing up to 50% of the flux 
see e.g.lWhite ct al. 1988; Church & Balucihska-Churcr, 
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2003al and references within). 

Van der Klis (1994) proposed that LMXBs switched 
from low to high accretion rate behaviou r at some con- 
stant f raction of the Eddington limit. IBarnard et al.l 
(|2003bh realised that if this were true, then BH LMXBs 
would be capable of exhibiting low-state variability and 
spectra at higher luminosities than neutron star LMXBs, 
as the Eddington limit is prop ortio nal to the mass of 
the accretor. IBarnard et al.l (|2004l ) found some empir- 
ical evidence for a transition at ~10% Edddington. We 
would therefore expect a LMXB containing a ~10 M Q 
BH to exhibit low-state characteristics at ~10 times the 
highest luminosity of low-states observed in 1.4 M Q NS 
L MXBs. 

iGladstone et al.l (|2007f) have systematically analysed 
observations of 40 disc-accreting neutron stars (atoll 
class LMXBs and millisecond pulsars) from the RXTE 
database, in order to follow their spectral evolution as 
a function of luminosity. They excluded sources with 
high line-of-sight or intrinsic absorption (equivalent to 
>3xl0 22 H atom cm" 2 ), as well as those sources which 
readil y exceed the Edding ton limit such as Sco X-l (sec 
e.g. IBarnard et al.l l2003af ). They found two different 
types of hard/soft transition: those where the spectrum 
softens at all energies, leading to a diagonal track on a 
colour-colour diagram, and those where only the higher 



energy spectrum softens, giving a vertical track. The di- 
agonal transitions were made at <10% Eddington, while 
vertical transitions occurred at <2% Eddington in the 
0.01-1000 keV band. 

Furthermore, GS 2023+338 (a.k.a. V4 04 Cygni) con- 
tains a BH primary with mass ~12 Mq (jShahbaz et al.l 
1994); it has exhibited low-state PDS and emission 
spectra at 2-37 ke V luminosities >3xl0 38 e rg s" 1 
([Miyamoto et al.l 11992b lOosterbroek et alJ 119971 ). This 
corresponds to a luminosity of >1.2±0.4xl0 38 erg s" 1 
in the XMM-Newton pass band (0.3-10 keV), for power 
law emission wi t h ph oton index 1.0-1.4, as found by 
iMivamoto et alJ (|1992f ). This translates as 8% Edding- 
ton for a 12 M Q black hole accretor, suggesting that spec- 
tral transitions in disc-accreting black holes may occur 
at similar fractions of the Eddington limit. 

We note that the M31 B H candidates identified by 
IBarnard et al.l (|2003bll2004f ) were contaminated by artifi- 
cial variability that was introduced when combining non- 
synchronised lightcurves f rom the different instr uments 
on board XMM-Newton ([Barnard et al.l l2007bf ). The 
current work employs corrected techniques described in 
IBarnard et all (|2007bD . 

3. KNOWN PROPERTIES OF BO 45, BO 135 AND BO 375 

3.1. X-ray properties 

The X-ray counterpart to Bo 45, has not been 
previously observed by XMM-Newton or Chandra. 
However, it was obse r ved b y Einstein in 1979-1980 
(|Trinchieri fc Fabbianol [1991 |) and by ROSAT in th e 
summers of 1991 a nd 1992 (jSupper et al.lll997l l200lT ). 
I Supper et al.l (|1997l ) designated the X-ray source as RX 
J0041. 7+4134, and found its 0.2-4.0 keV flux to have in- 
creased by a factor of —2 with respect to the Einstein 
observation, from 8.7+1.1 to 14.4+0.5 xl0~ 13 erg cm -2 
s — 1 ; this equates to a change in 0.2-4.0 keV luminos- 
ity from ~6x l0 37 to ~1.0xl0 38 erg s" 1 for a -780 
kpc distance (|Stanek fc Garnavichl ll998K For this es- 
timate, they assumed a 5 keV bremsstrahlung emission 
model, suffering absorption equivalent t o 7xl 20 H atom 
cm~ 2 , following Trinch ieri fc Fabbianol (|1991l ). However, 
they reported no variability between ROSAT observa- 
tions greater than 3a. We shall refer to this source as 
XBo45 for the remainder of the paper. 

The X-ray source associated with Bo 135 has been 
observed by both XMM-Newton and Chandra before; 
iTrudolvubov fc Priedhorskvl (|2004f ) have analysed these 
data as part of a spectral survey of M31 GCs. They 
found the 0.3-10 keV luminosity of Bo 135 to vary over 
3.3-4.1 xlO 38 erg s _1 , for power law spectra with photon 
indices varying over 1.48-1.66. We name this source 
XBol35. 

The X-ray source associated with Bo 375 has not been 
previously observed with XMM-Newton. However, it has 
been observed several times si nce 1979, with E i nstein , 
ROSAT, ASCA and Chandra; iDi Stefano et~aT1 (pOOJt ) 
have analysed many of these observations. They ob- 
tained good fits with absorbed power law models for the 
ASCA and ROSAT spectra, but not the Chandra ACIS- 
S spectrum; they obtained 0.5-2.4 keV luminosities of 
~ 1.4-2.3xl0 38 e r g s -1 , assuming a distance of 780 kpc. 
IDi Stefano et all (|2002f ) modeled the ACIS-S spectrum 
with a 0.80 keV blackbody and a power law with photon 
index, T = 1.67, for a 0.3-7 keV luminosity of 4.2xl0 38 
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erg s~ x . iTrudolvubov fc Priedhorskvl (|2004T ) re-analysed 
this Chandra observation, as well as an additional ACIS- 
S observation taken ^4 months later. They modeled the 
0.3-7 keV spectra with absorbed power laws, and found 
that the later observation was a factor ~2 brighter. We 
name this source XBo 375. 

3.2. Optical properties of the host clusters 

Bo 45, Bo 135, and Bo 375, the host clusters, have been 
classified as certain globular clusters in the " The Revised 
Bologna Catalogue (RBC) of M31 globular clusters and 
candidates" . They have been imaged with the HST and 
their natur e has been firmly c onfirmed (RBCV.3.4, Jan- 
uary 2008: lGalleti et al.ll2004L [2005. 20061 l2007f ). 

The integrated magnitudes and colours of the clusters 
can be found in RBC V.3.4 and can be used to provide 
estimates of the ages and masses of the clusters, in con- 
junction with reddening and distance and metal abun- 
dance values. Bo 45 has V = 15.78, V - I = 1.27, while 
Bo 135 has U = 16.04, V-I= 1.22, and Bo 375 has V = 
17.61, V — I = 1.02. We have adopted a distance mod- 
ulus o f (to - M) = 24.471 ± 0.0 35 ± 0.045 (784±13±17 
kpc, IStanek fc Garnavichlll998f k the first uncertainty 
is statistical and the second uncertainty is systemmatic. 
We use a reddening of E{B - V) = 0.10 ± 0.03, which 
is the aver a ge of the reddenings of all M31 clusters in 
iRich et all (|2005l ). Actually, for Bo 45, there are spe- 
ci fic reddening mea surements that have been tabulated 
in lRich et al. (2005), yielding an average of EiB — V) = 
0.10 ± 0.06, which is identical to the mean value we have 
adopted. Using these values, we have derived My = 
-9.0,-8.74,-7.17 and (V - I) = 1.13,1.08,0.88 for 
Bo 45, Bo 135 and Bo 375, respectively. The metal abun- 
dances of the three clusters are es timated to be i?=0.002 , 
0.0004 and 0. 0012 respectively dBellazzini et al.l ll995lh 
According to ISaraiedini et al.l (|2007l ) , the above values 
for (V - I) , M v and Z suggest that Bo 45 and Bo 135 
are at least as old as ^10 10 yr and at least as massive as 
~ 10 6 M Q , while the fainter optical luminosity and bluer 
color of Bo 375 allow it to be as young as 1.2 Gyr and to 
have a mass as low as 7xl0 4 M . 

In conclusion, both Bo 45 and Bo 135 seem to belong to 
the old massive globular cluster population of M31. This 
is not an unexpected result, since low mass X-ray binaries 
are found preferentially in luminous (massive) and red 
globular clusters (Kundu et al. 2008). However, Bo 375 
appears to be considerably less massive, and younger, 
than the other two GCs. 

4. OBSERVATIONS AND DATA ANALYSIS 

The 2006 December 26 and 2007 January 2 observa - 
tions of M31 with XMM-Newton (jJansen et all 120011) 
were conducted as part of the M31 extended survey pro- 
gram (PI W. Pietsch), which expands on the 2002 ma- 
jor axis surv ey to cover the e ntire optical D25 region 
of M31 (see IStiele et al1l2007lh In this work, we used 
the data from the European Photon I maging Cameras 
EPIC): EPIC-pn JStruder et al.ll200lh and EPIC-MOS 
Turner et al.ll2001lh We reconstructed the EPIC-pn and 
EPIC-MOS events files using the latest version of the 
XMM-Newton analysis package, SAS 7.1. A journal of 
observations for the three targets is presented in Table [TJ 

XMM-Newton observations are known to experience 
intervals of high background, so we screened these flares 



in the recommended manner. For the EPIC-pn, we cre- 
ated a 10-12 keV lightcurve with 100 s binning from the 
whole field of view, selecting only single events (PAT- 
TERN==0) from good pixels ( #XMMEA_EP); we re- 
jected all intervals with >0.4 count s _1 . Similarly, we 
created lightcurves for each EPIC-MOS detector, for all 
single events with energy >10 keV, with the equivalent 
EPIC-MOS filter (#XMMEA.EM). After checking that 
the flaring intervals were the same in the EPIC-MOS, we 
filtered the EPIC-pn and EPIC-MOS events files on the 
EPIC-pn flare filter; resultant good time exposures are 
provided in Table [T] 

For each source, we created a source extraction re- 
gion and corresponding background region that was at 
least as big as the source region, on the same CCD, at 
a similar off-axis angle and contained no point sources 
or extended emission. For the EPIC-pn events, we 
accepted PATTERN<=4, and FLAG==0, while for 
EPIC-MOS events we accepted PATTERN< = 12 and 
#XMMEA_EM. We created 0.3-10 keV source and back- 
ground spectra, along with the appropriate response files. 
We also created synchronised 0.3-10 keV lightcurves for 
the E PIC-pn, EPIC-MOS! and EPIC-MOS2 detectors 
(see iBarnard et al.l l2007al lbl). These lightcurves were 
then combined and their variability examined. 

To obtain astrometrically-corrected positions we se- 
lected the optical source s from the USNO- B1, 2MASS 
and Local Group Survey (jMassev et al.ll2006l ) catalogues, 
and checked that only one optical source was in the error 
circle of the corresponding X-ray source. We only ac- 
cepted sources correlating with globular clusters from the 
Revised Bologna Catalogue (V3.4) or with foreground 
star s, characterized by th eir optical to X-ray flux ra- 
tio dMaccacaro et al.lll988t ) and their hardness ratio (see 
IStiele et al.ll2008t ). For sources selected from the USNO- 
Bl catalogue, we used proper-motion corrected positions. 
We then used the SAS-task eposcorr to derive the offset 
of the X-ray aspect solution. 

5. RESULTS 

5.1. Lightcurve analysis 

The background regions of XBo 45, XBo 135 and 
XBo 375 had radii of 30", 50" and 40" respectively, while 
the source regions had radii of 20", 20" and 40"; these 
source regions were chosen to optimally balance encir- 
cled energy, background subtraction and source crowd- 
ing. We therefore normalised the background lightcurves 
by the background to source area ratios. The flare-free, 
normalised background lightcurves were stable at ^0.03 
count s _1 for all sources. 

We present the 0.3-10 keV lightcurves of XBo 45 and 
XBo 135 in Fig.[IJ binned to 400 s. XBo 45 is clearly vari- 
able; the r.m.s. variability is 9.2±0.9%, while the best fit 
line of constant intensity has a x 2 /dof = 356/108. Such 
variability is expected for low-state LMXBs, with most 
of the power at frequencies >1 Hz; hence, the true vari- 
ability of XBo 45 is expected to be higher than ~10%. 
We analysed the power density spectrum of XBo 45, in 
search of the characteristic broken power law PDS ob- 
served in low-state LMXBs. However, no significant 
power was detected in the PDS. This is likely to be due 
to faintness of the source. We note that XBo 135 exhib- 
ited no detectable variability; the best fit line of constant 
intensity yielded x 2 /dof = 115/108, and its r.m.s. vari- 
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ability was found to be 1.8±1.8%. The EPIC-pn+EPIC- 
MOS 0.3-10 keV lightcurve of XBo375 is shown in in 
Fig. [2 It is significantly less variable than XBo45 
(x 2 /dof = 153/86, 3.4±0.6% variability) despite having 
more than twice the intensity. 

5.2. Spectral analysis 

For each source we simultaneously fitted the EPIC- 
pn and EPIC-MOS spectra using XSPEC 11.3.21. Ta- 
ble Q] provides the number of net source counts in each 
detector from each source. For the models discussed 
below we assumed Solar abundances for our absorber 
(lAnders & Grevessd 11989ft. If we use the ISM abun- 
dances of lWilms et al.l ( 20001 ). and the TBabs absorption 
model in XSPEC, then we obtain column densities ~30% 
higher, but similar emission parameters and luminosities. 
However, the resulting fits are slightly worse. 

5.2.1. Fitting a power law emission model 

IWhite et al.l (1988) found the spectra of low luminosity 
(i.e. low state) LMXBs to be well described by inverse 
Compton scattering of cool photons on hot electrons in 
a corona, and a power law is a simplified representation 
of this model. Both XBo 45 and XBo 135 were well de- 
scribed by an absorbed, hard power law (see Tabled for 
parameters). The best fit power law models to the pn 
and MOS spectra of XBo 45 and XBo 135 are presented 
in Fig. [3] and Fig. [5] respectively. The quoted uncertain- 
ties correspond to 90% confidence limits. However, the 
best fit absorbed power law model to the XBo 375 must 
be rejected. 

We also modeled the spectra of XB o 45 and XBo 13 5 
with the COMPTT spectral model (see lTitarchuklll994l ). 
Our spectra could not sensibly constrain the tempera- 
ture, favouring ^50 keV. 

We calculated the luminosities of XBo 45 and XBo 135 
using the best fit power law emission models, and a dis- 
tance of 784 kpc; uncertainties in the distance will be 
discussed below. The unabsorbed 0.3-10 keV luminosi- 
ties of XBo 45 and XBo 135 were 2.46±0.09xl0 38 and 
4.76±0.11xl0 38 erg s _1 respectively. 

5.2.2. Adding a blackbody component 

A power law emisssion model is rejected for XBo 375; 
however, adding a second, blackbody, component yields 
a good fit, shown in Table[2j The pn and MOS spectra of 
XBo 375, along with this best fit model, are presented in 
Fig. The best fit, unabsorbed 0.3-10 keV luminosity 
for this source is 6.5±0.2 x 10 38 erg s _1 . The black- 
body component contributes 14±4%. We note that our 
values of YT and T are consistent with those found by 
iDi Stefano et all (|2002f) for the ACIS-S spectrum; how- 
ever, our preferred absorpt i on is a factor ~2 lower. We 
note that[Di Stefa no et al.l (2002) found the absorption 
to vary by a factor of ~3 between ROSAT, and Chandra 
observations, so this is likely to be real. 

For XBo 45, the fit is marginally improved by adding 
a blackbody component; however, the blackbody con- 
tribution to the flux could not be sensibly constrained. 
We infer from this that there is no significant blackbody 
contribution to the spectrum of XBo 45. 

Adding a blackbody component to the XBo 135 emis- 
sion spectrum improves the fit significantly. The best fit 



BB+PO model is shown in Table H The 0.3-10 keV 
luminosity for this model is 4.4±0.2x 10 38 erg s _1 ; the 
blackbody contributes 11±5%. 

The spectra of XBo 375 had twice as many counts as 
the equivalent spectra of XBo 45 or XBo 135. Hence we 
created an EPIC-pn spectrum for XBo 375 with a sim- 
ilar number of source photons to the EPIC-pn spectra 
of XBo 45 and XBo 135, to see if the blackbody compo- 
nent was still necessary. The best fit power law yielded 
X 2 /dof = 324/251 (g.f.p = 0.0013); hence, a blackbody 
component was still required; the best fit BB+PO model 
yielded x 2 /dof = 265/249 (g.f.p. = 0.22). These results 
show that the emission of XBo 375 is significantly differ- 
ent from that of XBo 45, where no blackbody is required. 

We then modeled the spectra of each source with ab- 
sorbed blackbody and disc blackbody models, as might 
be expected for a BH LMXB in the high state. The best 
fit blackbody models yielded x 2 /dof >4 for each source. 
The disc blackbody models fared better, but still must be 
rejected, with good fit probabilities ranging from 0.02 for 
XBo 135 to 10~ 7 for XBo 375. Hence, we may securely 
reject thermally dominated emission for each source. 

5.2.3. Fitting low-metallicity absorbers 

IDi Stefano et al.l (|2002f ) investigated the effect on mod- 
elling XBo 375 spectra of assuming a metallicity appro- 
priate to XBo 375 (6% Solar) for the absorber. Using 
the vphabs model in XSPEC, they set the H and He 
abunances to cosmic values, and fixed the abundances 
of higher elements to 6% Solar. Their modelling of the 
Chandra data was improved (x 2 /dof =1.29 for 206 dof, 
c.f. 1.53 for 209 dof for a standard absorber); however, 
neither fit was formally acceptable (good fit probability 
<0.003). Neither ROSAT nor ASCA spectra were im- 
proved by low metallicity fits. 

We therefore modeled our XMM-Newton spectra 
of XBo 45, XBo 135 and XBo 375, using metallicities 
of 0.11, 0.02 and 0. 06 respectively, as obtained by 
iBellazzini et al.l (|1995|) . We present the best fit models 
for each source in Table [3) We see that the fits are con- 
siderably worse in each case, with only the BB+PO fit to 
XBo 135 being a good fit. It is likely that the metallicitcs 
of the environs surrounding these sources were enhanced 
by the supernovae that produced the NS/BH accretors, 
and therefore the GC metallicities are inappropriate. 

5.3. The spectral evolution of XBo 4-5 

Of the three sources, only XBo 45 exhibited signifi- 
cant time variability. The observed intensity variation 
in XBo 45 could be due to one of two phenomena: varia- 
tion in the emission spectrum or variation in the absorb- 
ing material. To test these scenarios, we obtained two 
additional EPIC-pn spectra: one from a low intensity 
interval, and another from an interval of high intensity, 
represented in Fig. Q] by a solid line and a dashed line 
respectively. We applied the same response and back- 
ground files to these spectra as to the original XBo 45 
spectrum. The low intensity spectrum contained ^3000 
source counts, for a mean intensity of 0.424+0.008 count 
s -1 , while the high intensity spectrum contained ~4500 
source counts, with a mean intensity of 0.512+0.008 
count s . 

Simultaneously fitting the low and high intensity spec- 
tra with r and normalisation linked but free to vary, and 
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with the absorption (iVn) free to vary yielded an unac- 
ceptable best fit. We then linked JVh and T, varying 
only the normalisation; this produced an acceptable fit. 
Finally, we freely fitted each spectrum, and found that 
iVjj and r were consistent within 90% confidence lim- 
its for the two spectra. The best fits for each of these 
models are presented in Table |U It is clear that the 
variation is intrinsic to the X-ray source, rather than 
the abso rber; such variation is characteristic of low state 
LMXBs (|van der Klidfl99l . 

5.4. Are the X-ray sources associated with the GCs? 

Mass segregation is thought to concentrate binaries, 
which arc amongst the heaviest objects in globul ar clus- 
ters, to the centre of the cluster (e.g. iMevlan fc H cggic 
1997). The X-ray sources are each located ~1" from 
their respective GC positions, with positional errors of 
~1". Of the three GCs, only Bo 45 has sufficient q uality 
data to measure its structure. iBarmbv et al.l (|2007l ) used 
HST observations of Bo 45 to model its spatial charac- 
teristics, and found a core radius of 0.3" and a half-mass 
radius of 1". It is therefore impossible to confirm whether 
the X-ray sources are located in the central regions of the 
GCs. 

Following iMaccarone et all (|2007l ). we first calculated 
the probability that the X-ray sources have GC associ- 
ations purely through chance, and then calculated the 
probability that these sources are active galacic nuclei 
(AGN) that are coincident with the GCs. 

Both XBo45 and XBo 135 were observed in the NN1 
region (Table 1). Assuming a circular field of view with 
15' radius, we searched the literature for GCs and GC 
candidates within the NN1 observation, finding 16 GCs 
and 51 GC candidates. Hence the GC spatial density 
<3xl0 -5 per square arcsec. We found 30 X-ray sources 
with >100 photons in the EPIC-pn camera; therefore the 
probability of finding a GC within 1" of a bright X-ray 
sources is <9xl0 -4 , and the probability of two chance 
associations is <8xl0 -7 . 

XBo 375 was observed in the NS2 region (Table 1). 
We found 19 GCs and 39 GC candidates, giving a spa- 
tial density <2xl0~ 5 . There were 20 sources with >100 
EPIC-pn counts, and the probability of a GC coinciding 
with one of these by chance is <4xl0~ 4 . 

We estimated the probability t hat the X-ray source s 
were AGN by utilising the work of iMoretti et alJ (|2003h . 
who have modeled the AGN flux distributions in the 1-2 
keV and 2-10 keV bands. We found 0.33to;f AGN per 
square degree with 1-2 keV fluxes as high as the 1-2 keV 
flux of XBo 45 (4.4xl0~ 13 erg cm~ 2 s" 1 ), and 0.08±g;J| 
AGN per square degree with 2-10 keV fluxes as high for 
XBo 45 (2.3xl0 12 ergcm~ 2 s _1 ). Even if we assume that 
none of these 2-10 keV sources are included in the 1-2 
keV sources, and take the upper limits, then the spa- 
tial density of AGN as bright as XBo 45 is <7xl0~ 8 per 
square arcsec; the probability of such an AGN coincid- 
ing with a GC is <5xl0 -6 . Furthermore, XBo 135 and 
XBo 375 are more luminous than XBo 45, and they are 
even less likely to be AGN. 

We also note that any AGN, or indeed foreground star, 
would be very likely to distort the optical colours of the 
cluster. We therefore conclude that the X-ray sources 
are very probably associated with the GCs. 



6. DISCUSSION AND CONCLUSIONS 

We have examined the emission spectra and time vari- 
ability of three X-ray sources associated with GCs in 
M31, using the 2006 December 26 and 2007 January 2 
XMM-Newton observations. The emission of XBo 45 is 
well described by a pure power law with photon index 
~1.4, and is highly variable. This is consistent with a 
NS or BH LMXB in the low state (jvan der Klislll994h . 
but is not consistent wit h a BH LMXB in the high state 
or steep power law state (|McClintock fe Remillardir2006l 
and refe rences within), or a NS LMXB emitting at >10 38 
ergs " 1 (IWhite et alll 19881; IChurch fc Balucinska- Church! 
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Sta nek fe Garnavichl (1998) calculated a distance to 
M31 of 784 kpc, ±13 kpc of statistical error, ±17 kpc 
of systematic error. Combining these distance uncer- 
tainties adds further uncertainties in the source lumi- 
nosities of ±5%. Therefore, the 0.3-10 keV luminosity 
range for XBo45 is 2.5±0.2xl0 38 erg s" 1 , or 140±10% 
of the Eddington limit for a 1.4 M Q neutron star pri- 
mary. However, several LMXBs have been found to con- 
tain neutron stars with mass as high as —2.1 M (see e.g. 
iNiceet all [20051 lolel 120061) : XBo 45 has a 0.3-10 keV 
luminosity of —80% E ddington for such systems. Since 
iGladstone et al.1 (|2007h showed that transitions in neu- 
tron star LMXBs occur at <10% Eddington, we consider 
XBo 45 to exhibit low state behaviour at an apparent lu- 
minosity too high for a neutron star. We therefore iden- 
tify the accretor in XBo 45 as a BH candidate. We note 
that XBo 45 has been observed several times by the Ein- 
stein and ROSAT observatories over the last — 30 years, 
varying in luminosity only by a factor of —2. This be- 
haviour is consistent with that predicted for a GC BH 
binary formed by tidal capture (jKalogera et al.l l2004). 

In contrast to XBo 45, the observed two component 
emission of XBo 375 is consistent with a bright NS 
LMXB, but not a low state NS or BH LMXB, nor a BH 
in the high or steep power law states. Hence, we classify 
XBo 375 as a NS LMXB. Finally, the emission spectrum 
of XBo 135 is consistent with a pure power law with pho- 
ton index —1.6, but the fit is significantly improved by 
adding a blackbody component. Hence, XBo 135 is con- 
sistent with a NS or BH LMXB, and deeper observation 
is required for further classification. 

Of the thirteen bright X-ray sources in Galactic GCs, 
twelve have exhibited X-ra y bursts, confiming thei r na- 
tures as NS LMXBs (see IVerbunt fe Lewinl 120061 for 



a review). iPietsch fc Haberll ( 20051 ) conducted a survey 
of XMM-Newton observations of X-ray sources associ- 
ated with M31 GCs, looking for X-ray bursts. They 
found simultaneous pn and MOS detections of bursts in 
two sources, and several burst candidates that were de- 
tected in the pn only. Such bursts would be identifiable 
in the lightcurves of our target sources. However, X-ray 
bursts are t hought to be forbid den at luminosities > 50% 
Eddington (jLewin et al.lll993l ). hence we do not expect 
bursts from XBo 135 or XBo 375. 

XBo 45 exhibits low-state behaviour at luminosities 
-10 times higher than expected for 1.4 M NS LMXBs; 
XBo 135 may also. A composite of two or three bright X- 
ray sources would likely result in a different spectrum to 
that observed in XBo 45; hence, if XBo 45 were a com- 
posite, it would be more likely to be made from —10 
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NS LMXBs in the low state. iDi Stefano et all (|2002l ) 
calculated the probability for a GC to contain multiple 
bright X-ray sources: they argue that if p is the proba- 
bility of finding one bright X-ray source in a GC, then 
Poisson statistics dictate that probability for two bright 
X-ray sources should be p/2, while the probability for 3 
bright X-ray sou r ces sh ould be p 2 /6. For a p ~0.1-0.2, 
IDi Stefano et all (|2002l) predict 3-5 GCs with 2 X-ray 
sources and less than 1 with 3 X-ray sources in M31. 
It is therefore very unlikely that XBo45 combines the 
emission of ~10 low-state neutron star LMXBs. 

It is however possible that the emission from XBo45 
is anisotropic. It exhibits low-state spectra at lumi- 
nosities ^10-20 times higher than expected for neu- 
tron star LMXBs; hence it could simply be beamed 
by a factor of 10-20. In this case, the emission would 
be restricted to a small solid angle; one might expect 
to observe suc h beaming in ~5-10% of random ly ori- 
ented systems. iTrudolvubov fc Priedhorskyl (|2004f ) mod- 
eled the spectra of 43 M31 GCs. If we loosely class a 
low-state spectrum as a power law with F <1.7, then 
ITrudolvubov fc Priedhorskyl (J2004D found 20 GCs con- 
sistent with low-state spectra. Five of those, XBo5, 
XBo82, XBol35, XBol53 and XBo386 exhibited ap- 
parent luminosities >10 38 erg s -1 . If we include XBo45, 
then 6 out of 21 GC systems consistent with low-state 
spectra (including XBo45), have measured luminosities 
>10 38 , i.e. ^30%, a significantly larger fraction than 
expected from beaming. 

We note that the host cluster Bo 45, which contains 



a BH LMXB candidate, is significantly larger than the 
cluster Bo 37 5 (see Section 3.2 ) , whic h we think contains 
a NS LMXB. IDi Stefano et all (|2002h describe Bo 375 as 
not at all unusual, with parameters close to the median of 
M31 GCs. This suggests that Bo 45 (and also Bo 135) are 
particularly massive, and may therefore be more prone 
to forming BH LMXBs. Therefore we conclude that, 
unlike the Milky Way, at least one GC in M31 is likely 
to contain a black hole binary. 
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TABLE 1 

Journal of observations. For each object we give the position, XMM-Newton observation, good time exposure, and 

number of net source photons in the epic-pn, epic-mos1 and epic-mos2 spectra. 

Source X-ray Position Observation Date Exp pn cnt MOS1 cnt MOS2 cut 

XBo45 00 h 41 m 43.19 s +41°34'20.1" NN1 (00 h 41 m 52.8 s +41°36'36.0") 
XBol35 00 h 42 m 52.00 s +41°31'09.7" NN1 

XBo375 00 h 45 m 45.54 s +41°39'42.6" NS2 (00 h 45 m 50.4 s +41°30'44.7") 



2006-12-26 


40 ks 


16894 


7416 


6953 


2006-12-26 


40 ks 


15104 


5560 


6090 


2007-01-02 


32 ks 


31946 


15686 


13798 
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TABLE 2 

Best fit spectral models for. fitting 0.3-7.0 keV EPIC-pn and EPIC-MOS spectra from XBo45, XBo 135 and 

XBo375. The models are absorbed power law (PO), and absorbed blackbody+power law (BB+PO). For each model 

WE GIVE THE COLUMN DENSITY (nn/10 21 CM -2 ), BLACKBODY TEMPERATURE (kT/keV), PHOTON INDEX (I"), CONSTANT OF 

NORMALISATION FOR EPIC-MOS1 AND EPIC-MOS2 (N M 1 AND A r M2) ! X 2 /dof [GOOD FIT PROBABILITY] AND 0.3-10 KEV 

UNABSORBED FLUX. NUMBERS IN PARENTHESES INDICATE THE 90% UNCERTAINTY IN THE FINAL DIGITS. 



Source 


n H /10 21 cm- 2 


kT/keV 


r 


N M i 


A r M2 


X 2 /dof [gfp] 


-F0.3-10/IO 12 erg cm 2 


S" 1 


XBo 45 PO 


1.41(11) 




1.45(4) 


1.06(4) 


1.03(4) 


517/487 


0.17 


3.34(12) 




XBo 45 BB+PO 


1.46(19) 


1.23(19) 


1.57(6) 


1.06(3) 


1.03(2) 


501/485 


0.30 


3.2(2) 




XBo 135 PO 


2.76(12) 




1.56(3) 


1.12(3) 


1.06(3) 


467/435 


0.14 


6.45(15) 




XBo 135 BB+PO 


2.3(3) 


0.8(2) 


1.54(14) 


1.12(3) 


1.06(3) 


413/433 


0.75 


6.0(3) 




XBo 375 PO 


1.52(6) 




1.64(2) 


1.15(2) 


1.07(2) 


1180/1032 [9E-4] 9.17(17) 




XBo 375 BB+PO 


1.41(11) 


0.90(10) 


1.73(18) 


1.15(2) 


1.07(2) 


1110/1030 [0.19] 8.8(3) 
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TABLE 3 
Best fit models for fitting 0.3-7.0 keV EPIC-pn and EPIC-MOS spectra from XBo45, 
XBo 135 and XBo 375 using observed GC metallicities. For each source we give the 
column density (tih/IO 21 cm -2 ), temperature (kT/keV), photon index (r), constant of 
normalisation for EPIC-MOS1 and EPIC-MOS2 (N M1 and N M2 ), and x 2 /dof [good fit 
probability]. Numbers in parentheses indicate the 90% uncertainty in the final digits. 

Source n H /10 21 cm- 2 kT /keV T JV M i N M2 x 2 / dof [gfp] 

XBo45PO 0.188(12) ... 1.33(2) 1.06(2) 1.05(2) 567/487 [7E-3] 

XBo 135 PO 0.58(4) ... 1.39(2) 1.13(3) 1.06(3) 740/435 [3E-18] 

XBol35BB+PO 0.24(4) 0.69(2) 1.05(5) 1.13(3) 1.07(3) 436/433 [0.45] 

XBo375 BB+PO 0.131(14) 0.69(4) 1.42(3) 1.15(2) 1.07(2) 1110/1030[0.04] 
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TABLE 4 
Results from simultaneously fitting 0.3-7 keV EPIC-pn spectra from the high and low intensity intervals in the 

LIGHTCURVE OF XB0 45. FOR EACH MODEL, WE PRESENT ABSORPTION (JVh / 10 21 ), PHOTON INDEX (T), BEST FIT X 2 /dOF AND 
GOODNESS OF FIT. NUMBERS IN PARANTHESE REPRESENT 90% UNCERTAINTIES IN THE LAST DIGITS. 

Variable absorption Variable normalisation Free fitting 

Nn/10 21 cm- 2 T X 2 /dof [gfp] JVh/10 21 cm" 2 V X 2 / dof [g f P] Nn/10 21 V X 2 /dof [gfp] 

High 1.4(2) 1.51(8) 193/129 [1.9c-4] 1.5(2) 1.50(8) 140/129 [0.23] 1.5(3) 1.45(10) 125/127 [0.54] 

Low 1.5(3) 1.51(8) 1.5(2) 1.50(8) 1.4(4) 1.57(15) 
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Fig. 1.— Combined EPIC-pn+EPIC-MOS 0.3-10 keV lightcurves of XBo45 and XBol35, binned to 400 s. The lightcurves have been 
background-subtracted and screened for flaring. The lightcurve of XBo45 is clearly variable; its r.m.s. variability is 9.2±0.9%. However, 
XBo 135 shows no significant variability. The solid and dashed lines in the XBo 45 lightcurve indicate the intervals used for the low and 
high intensity spectra respectively, discussed in Section 5.3. 
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Fig. 2.— Combined EPIC-pn+EPIC-MOS lightcurve of XBo375 in the 0.3-10 keV band, binned to 400 s. Even though the count rate 
for XBo375 is 2-3 times higher than for XBo45 or XBol35, it is significantly less variable than XBo45. 
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Fig. 3. — Best fit power law model to simultaneously fitted 0.3-7 keV EPIC-pn (black) and EPIC-MOS (light and dark grey) spectra of 
XBo45. The upper panel shows the log-scaled spectra, while the lower panel shows the ratios of the observed to the expected flux for each 
energy bin 
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Fig. 4. — Best fit power law model to simultaneously fitted 0.3—7 keV EPIC-pn (black) and EPIC-MOS ( light and dark grey) spectra 
of XBol35. The upper panel shows the log-scaled spectra, while the lower panel shows ratios 



16 



Barnard et al. 

data and folded model 






CO 
-M 

c 

o 
o 



o 
o 



o 
o 



ro 
I 
O 



CM 

m 



ID 

o 




0.5 



1 2 

Energy (keV) 



Fig. 5. — Best fit power law model to simultaneously fitted 0.3—7 keV EPIC-pn (black) and EPIC-MOS ( light and dark grey) spectra 
of XBo375. The upper panel shows the log-scaled spectra, while the lower panel shows ratios 



